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ABSTRACT 

We present high resolution (R^100,000) L-band spectroscopy of 11 Herbig AeBe stars with cir- 
cumstcllar disks. The observations were obtained with the VLT/CRIRES to detect hot water and 
hydroxyl radical emission lines previously detected in disks around T Tauri stars. OH emission lines 
are detected towards 4 disks. The OH 2 n 3 / 2 P4.5 (1+,1-) doublet is spectrally resolved as well as the 
velocity profile of each component of the doublet. Its characteristic double-peak profile demonstrates 
that the gas is in Keplerian rotation and points to an emitting region extending out to ~ 15-30 AU. 
The OH emission correlates with disk geometry as it is mostly detected towards flaring disks. None 
of the Herbig stars analyzed here show evidence of hot water vapor at a sensitivity similar to that of 
the OH lines. The non-detection of hot water vapor emission indicates that the atmosphere of disks 
around Herbig AeBe stars are depleted of water molecules. Assuming LTE and optically thin emission 
we derive a lower limit to the OH/H2O column density ratio > 1 — 25 in contrast to T Tauri disks 
for which the column density ratio is 0.3 - 0.4. 

Subject headings: astrochemistry - molecular processes - protoplanetary disks - stars: pre-main se- 
quence, Herbig AeBe 



1. INTRODUCTION 

Young pre-main-sequence stars are surrounded by 
gas-rich dust disks that are the left over from the 
collapse of the molecular cloud core. Sub-micron 
sized dust grains grow to larger sizes wit h time (e.g. 
iBouwman et al.ll2001t iRodmann et al.ll2006h . This leads 
to the fo rmation of planetesimals and eventually plan- 
ets ( e .g. iWeidenschilling fc Cuzzil 119931 : iHenning et al.l 
12000: iBlum fc Wurml l2008f ). The evolution of the in- 
frared excess in pre-main-sequence stars (PMSs) suggests 
that most of the small dust grains in the disk disap- 
pear within 3-5 Myr from the collapse of the molecu- 
lar cloud and very little amount of dust is found be- 
ond 10 Myr (e.g. lHaisch etTal] 120011: IBouwman et al l 



2006t iHernandez et al l 120071 iRocc atagliata et al.l 120091: 
Pascucci fc Tachibanal I2010D . A similar timescale is 



found for the evolution of mas s accretion rates in 
disks (e.g. iMohantv et al.l 120051 : Ijavawardhan aTet al.1 
2fM ISicilijYA guilar etalT l20M llnglebv et al.l 120091 : 
Fedele et al.1 12010V The energetic radiation field of a 
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PMS star impinging onto the disk surface regulates the 
disk temperature, can ionize the gas, breaks molecular 
bonds and might lead to evaporation of the outer lay- 
ers of the disk. An important role for the formation of 
a planetary system as well as for the origin of life on 
Earth is played by water. Compared to other volatilcs 
(e.g. NH3, CO2, CH4) water has an higher condensa- 
tion temperature. Thus, within a protoplanetary disk, 
water is the first volatile to condense as temperature de- 
creases as a function of the radial distance from the star 
and vertical depth. In view of the cosmic abundance of 
hydrogen and oxygen, water is the most abundant ice. 
This sets a boundary (snow line) beyond which most of 
the molecular gas condenses onto ice. For a review on 
the snow line in the so l ar neb ula and in protoplanetary 
disks see, e.g., lPodolalj(|2010t ). Beyond the snow line the 
solid surface density, and per haps even the total surfac e 
density of the disk, increases ((Kennedy fc Kenvonl l2008') . 
This, in turn, speeds up the formation of gas giant plan- 
ets allowing them to form before the gas in the disk is 
dispersed. 

Water and hydroxyl (OH) emission have been de- 
tected in a number of protoplanetary di sks around 
young sun-like stars i n the near-infrared (|Carr et aT 



2001 iThi fc Biai2005l : ISalvk et al.l 12008: Mandcl l et al 



2008) as well a t mid- in frared wavelengths (10 - 40 /im) 
(ICarr fc Naiital 120081: ISalvk et al.l 120081 : iNaiita et al.l 
I2010t IPontoppidan et al.l I2010D and in the far-infrared 
with the Herschel Space Observatory ([Sturm et alj|2010t 
Ivan Kempen et aT]|201Cf ). Lines at different wavelengths 
trace different temperatures and hence different regions 
of the disk. While the near- to mid-infrared and high ex- 
citation far-infrared lines probe the molecular gas from 
the warm surface layers of the disk, the low excitation 
far-infrared lines are sensitive to colder gas (a few 100 K) 
locat ed closer to the disk midplane (e.g. iWoitke et al.l 
120091) . 

iMandell et al.l ([2008) performed high resolution 
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Table 1 

CRIRES observation log 



Star 


RA (J2000) 


DEC (J2000) 


Date observation 


Exposure (s) 


Airmass 


Calibrator 


UX Ori 


05:04:29.9 


-03:47:11.1 


2008-12-06; 


02:50:50 


1440 


1.2 


HD 36512 


HD 34282 


05:16:00.5 


-09:48:31.2 


2008-12-06; 


01:47:50 


1440 


1.1 


HD 40494 


CO Ori 


05:27:37.8 


+11:25:33.3 


2008-12-07; 


04:32:54 


720 


1.2 


HD 64503 


V380 Ori 


05:36:25.0 


-06:43:00.3 


2008 12 05; 


02:36:08 


720 


1.1 


HD 23466 


BF Ori 


05:37:12.9 


-06:35:07.5 


2008-12-07; 


05:15:18 


2400 


1.1 


HD 64503 


HD 250550 


06:01:59.5 


+16:30:50.7 


2008-12-06; 


05:34:49 


1080 


1.3 


HD 41753 


HD 45677 


06:28:17.3 


-13:03:18.2 


2008-12-06; 


07:31:57 


1080 


1.1 


HD 74195 


HD 259431 


06:33:04.6 


+10:19:16.6 


2008-12-06; 


06:22:57 


1080 


1.2 


HD 74280 


HD 76534 


08:55:08.6 


-43:28:01.1 


2008-12-05; 


07:04:53 


2400 


1.1 


HD 28873 


HD 85567 


09:50:28.2 


-60:58:02.9 


2008-12-06; 


08:22:59 


480 


1.3 


HD 98718 


HD 98922 


11:22:31.0 


-53:22:07.9 


2008-12-05; 


08:19:43 


720 


1.3 


HD 39764 



(R ^27,000) L-band spectroscopy towards two Herbig 
AeBe stars (AB Aur and MWC 758). They did not find 
evidence of hot water vap or emission in these two disks. 
iPontoppidan et all (|2010ft searched for colder molecular 
emission in the mid-infrared the Spitzer Space Telescope. 
In contrast to T Tauri stars, the mid-infrared spectrum 
of Herbig AeBe stars is poor in molecular emission and 
they hnd only tentative H2O and OH emission lines 
towards some Herbig AeBe stars. 

The different molecular emission between T Tauri and 
Herbig AeBe stars might be due to photochemistry which 
controls the excitation and dissociation of molecules. 
Aiming to test this hypothesis, we performed a deep 
search for hot water and hydroxyl radical vapor emis- 
sion in Herbig AeBe stars with circumstellar disks. We 
performed ultra high resolution (R ^100,000, Av ~ 3 km 
s^ 1 ) L-band spectroscopy of 11 Herbig AeBe stars with 
CRIRES at the VLT. This spectral resolution allows us 
to resolve the velocity profile of the molecular emission 
lines and thus determine the radii over which the emis- 
sion arises. 

Observations and data reduction procedures are pre- 
sented in Sec. [2j In Sec. [3] we describe the immediate 
results of the survey: the detection of the OH 2 n 3 / 2 P4.5 
(1+,1-) rovibrational line and the non-detection of water 
rovibrational lines. In Sec. 0] we present the analysis of 
the OH lines and discuss the implications of our findings 
in Sec. [5] Finally, we summarize our results in Sec. H3 

2. OBSERVATIONS AND DATA REDUCTION 

High resolution L-band spectroscopy was obtained on 
the nights of 4, 5 and 6 December 2008 with the ESO's 
VLT cryogenic h igh-resolution infra red echelle spectro- 
graph (CRIRES, iKaeufl et al.1l200l at the Paranal ob- 
servatory in Chile. A slit width of 0'.'2 was adopted. In 
order to center the target inside the narrow slit we used 
the adaptive optics system and the targets themselves as 
reference stars. This procedure reduces slit losses result- 
ing in a higher signal-to-noise ratio and higher spectral 
resolution. Unresolved sky emission lines can be used to 
determine the actual spectral resolution. The full-width- 
half-maximum (FWHM) of the OH sky lines measured 
in the raw frames is ~ 3 km s _1 (or R ^100,000) as ex- 
pected from the nominal resolution. The spectra were 
recorded by nodding the telescope along the direction 
perpendicular to the slit (oriented along the parallactic 
angle) with a nodding throw of 10". The spectra cover 
the wavelength range between 2861 — 2936 nm (order 19, 



X re f = 2909.6 nm) thus covering sev eral ro-vibrational 
lines of H 2 Q detected in come ts (e.g. Idello Russo et al.l 
120041: IPello Russo et al l 120051 ) and in protoplanetary 
disks (|Salvk et al.1 12001 as well as ro-vibrational lines 
of OH detected i n disks by iMandell et "all (|2008D and 
ISalvk et al.1 (120011 . 

The data were reduced using the ESO CRIRES pipeline 
v.ll.c£3 following a standard approach: first the 
CRIRES frames are corrected for flat field, dark and bad 
pixels; second the frames at a given nodding position 
are combined together and the two master frames (one 
for each nodding position) are then combined together 
after correcting for the spatial offsets (due to the nod- 
ding procedure); finally, the spectrum is extracted using 
a rectangular mask. The sky emission lines are used to 
determine the wavelength dispersion solution. To prop- 
erly correct for atmospheric telluric absorptions the spec- 
tra were bracketed with observations of standard stars of 
early spectral type (see Table [1]). These spectra are di - 
vided by the stellar atmosphere models of lKuruczl (|1979f ) 
to determine the instrument response function and at- 
mospheric transmission curve. The optical depth of the 
telluric lines is scaled to the depth of the science tar- 
get. Finally, each science spectrum is divided by the re- 
sponse function. The spectral regions heavily absorbed 
(atmospheric transmission < 0.5) are not used in the rest 
of the analysis. For the flux calibration we scaled the 
CRIRES s pectra to the observed L-band magnitude mea- 
sured by Ide Winter etHI pool V380 Ori, HD 76534, 
UX Ori, HP 25055 0, HD 259431, HD 45677, BF Ori), 
iMalfait et ail (Il99l HP 9892 2, HD 34282, HD 85567) 
and lMorel k Magnenad (| 19781 . CO Ori). 

Given the high temperatures of the stellar photosphere 
of Herbig AeBe stars, the L-band spectrum is relatively 
free of stellar photospheric absorption lines. The spectra 
of all the stars observed with CRIRES are continuum 
dominated and have high signal-to- noise ratio (S/N > 
100). Such an high ratio allows to detect weak emission 
lines down to a few percent of the continuum flux level. 
The exposure times vary from object to object and range 
from ~ 10 min (for targets of L ~ 4 — 6 mag) to 40 min (L 
~ 7 mag) . The log of the observations is given in Table 

m 

3. IMMEDIATE RESULTS 

11 http://www.eso.org/sci/data- 
processing/ software / pipelines / index, html 
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Table 2 

Properties of the program stars 



Star 


Sp.Type 


Age 


A v 


D 


Luv 


M acc 


Group 






[Myr] 


[mag] 


[PC] 


[Lq] 


[1O" 7 M yr- 1 ] 




UX Ori 


A3 


4.5 a 


0.8 C 


460 g 


2.4 


0.66 1 


II 


HD 34282 


AO 


7 b 


0.6 C 


160 e 




<0.2J 


1 


CO Ori 


P7 


< 0.12 c 


2.2 C 


460 g 


2.9** 


0.9 k 


II 


V380 Ori 


Al 


<0.01 c 


1.43° 


430 g 


18.8 


25 1 


1 


BF Ori 


A5 


3.2 a 


0.37 f -0.7 c 


460 g 


0.6-1.3 


0.87 1 


II 


HD 250550 


B7 


0.25 c 


1.17 c 


700 g 


160 


0.16 1 


1 


HD 45677 


B2 






>500 h 


>65 




11 


HD 259431 


B5 


<0.01 c 


0.63 f -1.62 c 


800 g 


225-2200 


7.8 1 


1 


HD 76534 


B2 


>0.5 d 


0.80 f 


870 g 


1020*** 




1 


HD 85567 


B5 


<0.01 c 


2.23 c 


1500 c 






II 


HD 98922 


B9 


<0.01 c 


0.54 c 


1000 c 


545 


17.4J 


11 



* Integrated between 1100 - 2430 A, (IUE spectrum, Valenti et al. 2003) and multiplied by a factor 10 - 4 ' A i"°. 
Where A1770 is the extinction at A = 1770 A measured from Av using the R = 3.1 extinction relation of Cardclli et al. 
119891) . 

** Missing short wavelength IUE spectrum, integrated between 1850 - 2430 A. 

Missing long waveleng th IUE spectrum, integrated between 1150 - 2000 A. 
a iMont csinos et aT| 1120091) : 
b IMerfn et alj <2004l) : 
c IManoi etahl 12006T ); 
d IMarti n-Zaidi et al. 1 2 008T); 

van den An ckcr et al. lfl99Sl ); 
f I Valenti etaLl 12003]'); 
g IHillenbrand et al.l ljl992T) : 
h Ide Winter fc van den Anckerl 1 119971) ; 
1 Donehew et al. submitted; 
j IGarcia Lopez et al.l 12006]) ; 
k ICalvet et~aTl p004D : 



We have searched for H2O and OH rovibrational lines 
previously detected in protoplanetary disks around sun- 
like stars, comets, as well as transitions that are not yet 
detected but reported in recent synthetic molecular line 
lists. In the following section we summarize the main 
findings for H2O and OH. 

3.1. Detection of OH lines 

The rovibrational OH 2 n 3/2 P4.5 (1+,1-) at 2.9345^ 
transition is detected towards 4 disks (V380 Ori, HD 
250550, HD 259431 and HD 85567, Fig. H}. The high 
spectral resolution of CRIRES allows us to resolve not 
only the doublet but also the velocity profile of each tran- 
sition. Modeling of the line profiles and estimates of the 
radii traced by this emission are presented in Sec. El 
The profile of the 2 n 3 / 2 P4.5 (1+,T) line varies from 
object to object (Figs. H©. In the case of V380 Ori 
each line shows a characteristic double-peak profile as 
one would expect if the gas is in Keplerian rotation. In 
the case of HD 250550 the lines are narrower reflecting 
the lower inclination of the disk (Table|3]). Finally, the 
OH line in HD 259431 and HD 85567 is very broad and 
the two transitions of the doublet are blended even at the 
high resolution of CRIRES. In the case of HD 259431, 
multiple peaks are visible hinting to multiple velocity 
components. Two peaks are detected (above 3a) at the 
sides of the line at 2.9340 /im and 2.9352 fim respectively. 
The position of these peaks is symmetric to the center of 
the emission and correspond to a velocity of ± 47 km s _1 . 
The OH emission is very broad (FWHM = 90 km s _1 ) 
and the sides of the emission might be contaminated by 
the strong telluric absorptions at 2.9335 /jm and 2.9355. 



We used different telluric standard stars (taken at differ- 
ent airmass) to check the reliability of the line profile. In 
all cases we found the two high velocity wings above the 
3cr level. Further confirmation of the presence of these 
wings comes from a posteriori analysis of the line profile 
(Sec. I4.2p : if the OH emission is made of a single veloc- 
ity component, the central part of the emission would be 
stronger than the sides due to the blending of the line 
(see e.g. dotted-dashed line in Fig. [5]). The analysis of 
the line profile confirms the presence of multiple velocity 
components and the reliability of the high velocity peaks. 
At higher velocity, where the atmospheric transmission 
function drops below 70%, the emission is heavily cor- 
rupted by telluric absoprtion hence we excluded these 
regions from the analysis of the line. 
Within the spectral range covered by our observations 
there are other two transitions of the OH 2 n 3 / 2 P branch 
(J = 2.5, 3,5) as well as two transitions of the OH 2 Hi/2 P 
branch (J = 2.5, 3.5). Unfortunately all four transitions 
fall either in highly saturated telluric absorption regions 
or inside the inter-chips gaps hence we cannot compute 
stringent line flux upper limits. 

3.1.1. OH column density 

The energy level of the ro-vibrational transitions of 
OH might be populated either thermally or radiatively 
(fluorescence). In the first case the collisions with atoms 
and molecules are the dominant excitation mechanism 
and the emission is only dependent on the temperature 
and density of the gas as we will show below (LTE case). 
In the case of fluorescent emission the ro-vibrational lev- 
els are populated by absorption of near-infrared photons 
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Figure 1. Portion of the CRIRES spectra for V380 Ori, HD 250550, HD 259431 and HD 85567. In each panel the upper line (red) is the 
extracted target spectrum along with the atmospheric transmission curve (blue line). The lower plot of each frame is the spectrum cleaned 
by the telluric absorption lines multiplied by a factor of 10. The dashed line shows the 3a detection level. For each star the OH P4.5 
(1+,1-) doublet is detected and spectrally resolved. The vertical lines show the positions of some of the strongest H2O rovibrational lines. 
The dotted line indicates the expected line height of the H2O (001-000) [1166 - 1267] transition in the case of OH/H2O column density 
ratio equal to unity. 
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Figure 2. OH P4.5 doublet towards V380 Ori. The red line is 
the best fit model (Sec. 14.21 Tabic [3]l assuming Keplerian rotation. 
The asymmetry of the line profile indicates a deviation from the 
Keplerian motion of the gas (see Appendix) . 



(in the electronic ground state) or absorption of UV pho- 
tons which excite the electronic states with following 
cascade in the electronic ground state. UV fluorescent 
pumping in proto planetary disks has been in vestigated 
for the first time bv lBrittain et al.l ([2003, 2007) to explain 




Figure 3. \ contours in the a — R parameter space for V380 
Ori assuming i = 30° and Ri = 2AU. The best fit value is shown 
as a black dot. 



the CO emission from the disk around the HAeB e stars 
HD 141569 A and AB Aur. IMandell et all ([20081 ) imple- 
mented UV and IR fluorescent pumping for the OH ro- 
vibrational lines detected in the disk around AB Aur and 
MWC 758. They find that a fluorescent model can ex- 
plain the observed OH line intensities as well as OH ther- 
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Figure 4. Same as Figure[2]for HD 250550 
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Figure 5. Same as Fig. Q] for HD 259431. The dashed and dot- 
dashed lines are the the optically thick and disk's atmosphere com- 
ponents respectively. 

mal emission. The difference between thermal and fluo- 
rescent emission is in the inferred OH column densities, 
which are higher in the LTE case. Since our observa- 
tions do not allow us to discriminate between LTE and 
fluorescent emission, we assume in the following analy- 
sis that the OH gas is in LTE, the emission is optically 
thin, and compute the total number of OH molecules 
from the measured line flux. For the g as temperature 
we tak e a value of 700 K as estimated by iMandell et al.1 
(2008) for the Herbig AeBe stars AB Aur and MWC 
758. With thes e assumptions the OH line intensity is 
(lHerzberglll9"50h 

I(OH) =N n (OH)ln>cA nm (1) 
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Figure 6. \ 2 contours in the eti — R,2,o parameter space for HD 
259431 assuming i = 50°, ai = 0, Ri,j = 0.9 AU and Ri, = 
R2,i = 1.4 AU. The best fit value is shown as a black dot. 
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Figure 7. OH P4.5 doublet towards HD 85567. 

with I(OH) the OH line luminosity (= 4 n d 2 Fqh) with 
d being the distance to the star and Foh the measured 
OH flux), N n the number of molecules in the upper state 
n, h Planck constant, c the light velocity, v the wavenum- 
ber of the transition (cm -1 ), and A nm the Einstein coef- 
ficient of the 2 n 3 / 2 P4.5 transition. In LTE the quantity 
N n relates to the total number of OH molecules (N) as 

Nn(OH) = gj exp(-E n /kT) 
N(OH) Q(T) 1 ' 

where gj ( = 2J+1) is the degeneracy of the state, E n the 
energy of the upper state, k the Boltzmann constant and 
Q(T) the partition sum. Substituting N n from Eq. [2] in 
Eq. Q] we have 
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N(OH) 



I(OH)Q(T) exp(E n /kT) 
hc£>A nm (2J + 1) 



(3) 



Table 3 

Stellar and disk parameters 



The transition para meters and the partitio n sum are 
taken from HITRAN (jRothman et all 12009ft . The to- 
tal number of molecules can be converted into a vertical 
column density (N co i) by dividing it for the area of the 
OH emitting region, 7r(R 2 - R?), with R; and R being 
inner and outer radius of the OH emitting region (see 
Sec. |53). When the OH P4.5 doublet is not detected 
and the emitting region is not known we assume Ri = 
0.4 AU and R Q = 10 AU and a line width of 30 km s" 1 
(based on our results of the OH line detections). The 
values of N and N co i are listed in Table [4] We note that 
increasing the temperature to T(OH) = 1000 K decreases 
N (hence N co i) almost by an order of magnitude. 

3.2. Non-detection 0/H2O lines 

In order to search for faint water vapor emission we 
fir st created a l i ne lis t using the water line list catalog 
of iBarber et al.l (|2006ft . The latter is the most complete 
and accurate water line list in the literature. We used 
their BT2 code to produce the synthetic spectrum and 
adopted the same temperature as that for the OH gas 
(700 K). Some of the strongest water lines detected by 
iSalvk et al.l ((2008ft in the T Tauri stars AS 205 A and 
DR Tau lie around 2930 nm close to the OH 2 n 3/2 P4.5 
doublet we have detected in 4 disks and hence in a clean 
part of the spectrum. For comparison Fig. [T] shows 
the positions of th e strongest water emission lines from 
ISalvk et all (|2008ft . None of the 11 Her big AeBe stars ob- 
served with CRIRES show evidence of hot water vapor 
emission in the L-band spectra. 

We compute the 3-cr upper limit of the line flux for the 
line H 2 (001-000) [11 66 - 12 67 ] at 2.931 fim detected 
in T Tauri stars. The upper limit is measured as the 
product 3 x a x AA, where a is the standard deviation 
of the spectrum and AA is a charatcristic line width of 
30 km s _1 or 2.9 x 10~ 4 ^im (of the same order of the OH 
line width). The line flux upper limits for the program 
stars are given in Table |4] Typical upper limits are of 
the order of 10~ 16 — 10~ 15 erg cm ~ 2 s _1 . The water lines 
detected by ISalvk etail (1200* Fl have fluxes of the order 
of 10" 15 - 10" 14 ' 



' erg cm 2 s 1 . 



3.2.1. H2O column density 

With the same assumptions of Sec. 13. 1 . T1 (optically thin 
emission, LTE, T = 700 K) we can estimate the upper 
limit to the total number of molecules and column den- 
sity of H2O. In the case of water, Eq.|3]must be modified 
to take into account an additional degeneracy associated 
with the nuclear spin (gi = 3 for J odd and gi = 1 for 
J even). We compute the upper limits of N(H 2 0) and 
N co i(H 2 0) using the line flux upper limits of the (001- 
000) [1166 _ 1267] transition. The results are listed in 
Table |4] We note that also in this case the number of 
molecules decreases by almost an order of magnitude by 
increasing the temperature to T(H 2 0) = 1000 K. 

12 Estimated from the Fig. 2 of ISalvk et"all pOOSft 
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In the case of HD 259431 the model is the sum of two velocity 
components. The subscript "2" refers to the second velocity 
c omponent 

a IHubrie etaLI 12009T) : 
b IHernandez et al.l {2004ft ; 
c IKraus et al.l l|2008l i: 

d Fixed; 

e Constrained to be > Rj, o 

3.3. OH/H 2 line ratio 

From the com parison of Fig. [T] with Fig. 2 of 
ISalvk et al.l (|2008ft it is evident that the OH/H 2 line ra- 
tio is higher in Herbig AeBe disks than in T Tauri disks. 
For the Herbig stars with detected OH emission lines, 
we have that the OH/H 2 line flux ratio is > 1.5 - 3(0- 
This translates into a lower limit to the OH/H 2 column 
density ratio of 1 - 25 (Table|4j. For comparison we show 
in Fig. □ the expected line height of the H 2 (001-000) 
[1166 ~ 1267] transition in the case of OH/H 2 column 
density ratio equal to unity (dotted line). The inten- 
sity of the OH P4.5 line for the four Herbig AeBe stars 
(Table g) is of the order of 10~ 15 - 10~ 14 erg cm" 2 s -1 , 
simila r to the OH and H 2 Q line flux measured in T Tauri 
disks (ISalvk et al.ll2008l) . However, in the case of T Tauri 
disks the OH/H 2 line flux ratio is ~ 1 and the column 
density ratio is 0.3 - 0.4, that is, the OH/H 2 column 
density ratio is ~ 3—25 times larger in Herbig AeBe disks 
than in T Tauri disks. This suggest that, in contrast to 
T Tauri disks, water vapor is less abundant than OH in 
the disk atmosphere of Herbig AeBe stars. 

4. ANALYSIS 

In order to understand the pattern of detections and 
non-detections we investigate whether there is any cor- 
relation between stellar and disk properties and the sur- 
veyed emission lines (Sec. 14. ip . We also model the pro- 
files of the resolved OH lines and provide estimates for 
the radial distances traced by the L-band OH emission 
in disks around Herbig AeBe stars (Sec. 14. 2\ . 

4.1. OH lines and disk properties 

A noticeable distinction appears when detections are 
compared with the disk geometry: the majority of the 
sources with detection of OH emissi o n ap pears to have 
a flared disk geometry. M eeus et al.l (|200lD have shown 
that the ratio of the far-infrared to the near-infrared flux 
is sensitive to the geometry of the disk. Far- infrared 
bright stars belonging to the group I of iMeeus et al.l 

13 We refer here to the total flux of the OH P4.5 doublet, i.e. 
twice the value listed in the fourth column of Table [4] 
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(2001) are thought to have a flared disk geometry, that 
is the disk scale height (H) and opening angle (H/R) 
of the disk increases with the distance from the star 
(R). Group II sources are thought to have a flat, self- 
shadowed geometry, due to the disk inner rim casting 
a shado w at larger disk r adii (iDullemond fc Dominikl 
120041) . Ivan Boekel et all (|2003j) found that the two 
groups occupy two different regions in the IRAS myi — 
TO60 color versus £nir/£ir diagram, where Lnir is the 
integrated luminosity as measured by the J, H, K, L and 
M photometry and Lir is the integrated luminosity mea- 
sured by the IRAS 12 25 and 60 ^m fluxes. Following 
Ivan Boekel et all ([2003) we provide in Table [5] the group 
of the program stars. Interestingly, three out of four stars 
for which OH emission is detected (V380 Ori, HD 250550, 
HD 259431) belong to group I while only one source 
showing OH emission (HD 85567) belongs to group II. 
We note that the other two HAeBe stars with detected 
OH emission in the literature (AU Aur and MWC 758 
iMandell et al.l I2008D also belong to gr oup I. We discuss 
the origin of such a correlation in Sec. 15.51 

4.2. OH line profiles 

In this section we present an analysis of the OH P4.5 
(1+,1-) line profiles for V380 Ori, HD 250550, HD 259431 
and HD 85567. Figs. [2][5]show the velocity profile of the 
doublet for the 4 stars from which we can estimate the 
radial distribution of the molecular gas. We assume that 
the line intensity follows a power-law distributio n as a 
funct i on of the radial distance from the star (e.g. iSmakl 
fl98lt iCarmona et al.l [20071: Ivan der Plas et al.l l2009fi ""of 
the form 

Ioh(R) = Ioh(Ro) • (R/R r Q (4) 

where R is the distance from the star and Ioh(Ro) is the 
intensity at the inner radius. We also assume that the 
OH line is optically thin. If the line is optically thick 
we should include i n the analysis the cont inuum and line 
optical depths (e.g. lHorne fc Mar sh 1986). Here, we are 
interested in determining the inner and outer radius of 
the OH emitting region and a complete analysis of the 
line emission properties is beyond the scope of this paper. 
The radial profile of the OH line is converted into a veloc- 
ity profile assuming that the gas is in Keplerian rotation. 
In addition the model line is convolved with a velocity 
width v = \/vm 2 + vth 2 where v; n is the instrumental 
broadening (~ 3 km s" 1 ) and Vth (= -\/2kT/moH) is the 
thermal broadening of the line (where k is the Boltzmann 
constant, T is the gas temperature and hioh is the mass 
of OH). We assume a temperature of 700 K as found by 
IMandell et al.l (|2008| ) which corresponds to vth ~ 0.8 km 
s _1 . For the Keplerian rotation and line convolution we 
us ed the IDL c odes "keprot" and "convolve"___| described 
in lAcke et all (|2005[ ). 

With the assumption of Keplerian rotation and of the 
radial profile in Eq. [4j the velocity profile of the emis- 
sion line is a function of 4 parameters: the power law 
intensity exponent (a) , the disk inclination (i) , the inner 
and outer radii of the OH emitting region (Ri,R ). The 
Keplerian velocity is Vk cp = \J GM„,/R, hence the inner 

14 http:www.ster.kuleuven.ac.be~bramdark_themework.html 



radius determines the maximum velocity in the line pro- 
file. Because the disk inclination is not known for our 
sources (apart from HD 259431) we fitted the OH P4.5 
doublet by minimizing the reduced \ 2 (x 2 ) between the 
observed and model velocity profile 

_ 9 1 v-^ /model — observed \ 2 , , 



with d degree of freedom. In our fitting procedure we 
leave free only three parameters: i, R; and R and repeat 
the model fitting for five different values of a (0, 1.0, 1.5, 
2, 2.5). The best fit parameters are listed in Table [3] 

4.2.1. Determination of the outer disk (R a ) 

The measure of R requires further explanation. The 
contribution from the outer disk to the line intensity is 
reduced because of the power-law distribution (Eqj4]). 
Moreover, the fit of the line profile might be degenerate 
if the parameters are not independent. As we saw above, 
Ri is uniquely determined by the maximum velocity of 
the line. In a Keplerian motion, the disk inclination and 
outer radius are responsible of the characteristic double- 
peak profile: for a given ring of gas at a distance R from 
the central star the two peaks of the line are separated by 
the quantity Vk op ■ sin(i) cx R -1 / 2 • sin(i). If the disk in- 
clination is known, it is possible to unambiguously deter- 
mine R . The power-law index (a) controls how steeply 
the intensity intensity decreases with radius (hence veloc- 
ity) and affects only slightly the p eak-to-peak separation 
as shown in Fig. 1 of Smak (198lj). To check the robust- 
ness of the fit we first determine the best fit parameters, 
then we fix the disk inclination and inner radius and 
compute the % 2 surface on the a — R parameter space. 

Notes on individual targets: 

V380 Ori. The OH emission extends from — 2 AU to ~ 
15 from the central star (Fig. ^ and we estimate a disk 
inclination of ~ 30°. Fig. [3] shows the \ 2 of the fit in the 
a — R plane. The single peaked x 2 distribution shows 
that the determination of R is robust. The OH P4.5 line 
profile is asymmetric with the blue peak stronger than 
the red peak. The model in Fig. [5] has been manually 
modified to account for the excess flux in the blue-shifted 
peak by multiplying the latter by a factor of 1.2. The 
origin of the asymmetric line profile is discussed in the 
Appendix. 

HD 250550. In this case the line is top-flat and the two 
peaks are not visible. This is likely due to the low incli- 
nation of the disk in the plane of the sky (disk almost 
face-on). In this situation it is not possible to determine 
unambiguously R so we fix it to a value of 100 AU. The 
best fit is found for an inner radius of ~ 0.4 AU and an 
inclination of 10°. 

HD 259431. The OH line profile in this case is more com- 
plicated and appears to have multiple components. We 
note that we were not able to fit the OH line profile with 
a model made of a single velocity component as in the 
previous cases. Multi- wavelength interferometric obser- 
vations reveal the presence o f optically thick g as within 
the dust sublimation radius (jKraus et al.l 120081 ). In par- 
ticular, the near-infrared continuum is dominated by op- 
tically thick gas that is accreting onto the star, while the 
mid-infrared continuum arises from the passively irradi- 
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ated disk atmosphere at lar ger radii. Similarly, the anal- 
ysis of the H 2 FUV lines (jBouret et al.ll2003f) suggests 
that the molecular hydrogen spectrum (seen in absorp- 
tion) has multi-temperature components. In particular 
they find that an hot component (T ~ 1300 K) comes 
from optically thick gas in the vicinity (> 0.5 AU) of the 
star. This is in good agreement with the detection of the 
high velocity wings in the OH line profile (Sec. [3]). To 
account for the presence of optically thick gas inside the 
dust truncation radius we modify the standard model in- 
cluding a second component (see Fig. [S]). Thus the free 
parameters of the fit are now 5: disk inclination, inner 
and outer radius of the optically thick component (Ri,;, 
Ri,o) ; inner and outer radius of the passively irradiated 
disk component (R2,i, R 2 ,o)- We adopt a disk in clination 
of 50° based on the results of IKraus et all ([20081 ) . For the 
thermal broadening of the line we assume Ti=1300K (= 
T(H 2 )) and T 2 = 700 K (from IMandell et al.ll200ll for 
the high and low velocity component respectively. Fi- 
nally, we assume a constant radial dependence of the 
line intensity from the radius (a ~ in Eq. [4j for the 
high velocity component. Our best model parameters 
are in good agreements with the results of the interfero- 
metric observations: the first component extends from ~ 
0.8 AU to ~ 1.4 AU and the second components extends 
from ~ 1.4 AU to — 25 AU from the star. We note that 
our fitting routine tends to produce a low value of R 2 .i 
(lower than Ri. ), thus we constrain R 2 ; to be > Ri, - 
Fig. [5] shows the x 2 of the fit in the a — R 2 plane (after 
fixing all the other parameters). The value of R 2i0 is not 
degenera te with a. 

Recently iBagnoli et al.l (|2010f) found a similar result for 
the [O I] line at 6300 A. High spectral resolution obser- 
vations revealed the presence of an high velocity compo- 
nent between ~ 1 - 2 AU and a low velocity component 
peaking around 20 AU. Compared to the [O I] line profile, 
the OH high velocity component is clearly double-peaked 
and allows us to constrain better the outer radius of the 
high velocity component (Ri, )- Compared to the [O I] 
6300 A line, the OH emission detected here likely arises 
from deeper layers (higher Ay) in the disk. 
HD 85567. The OH emission is broad (of the order of 
~ 100 km s -1 ) and suggests the presence of high veloc- 
ity molecular gas. As in the case of HD 259431, the 
high velocity OH emission might originate in an opti- 
cally thick gas inside the dust sublimation radius. There 
are no estimates of the disk inclination in the literature 
and given the low signal-to-noise of the spectrum it is 
hard to estimate the radial extent of the OH emitting 
region. Further observations are needed to measure the 
properties of the OH emission from this disk. 

5. DISCUSSION 

The main findings of Sec. [3] & Sec. 0] are: i) the non- 
detection of water vapor emission lines, ii) the detec- 
tion of OH emission lines from the disk atmosphere of 
Herbig AeBe stars and iii) the correlation between OH 
emission and disk flaring. In contrast, observations of 
T Tauri stars have shown that there exist a hot (500- 
1000 K) water layer in their disk atmosphere with den- 
sities sim ilar to the OH layer we detect in Herb ig AeBe 
stars (e.g lCarr fc Naiitall2008HSalvk et alj|2051 . In this 
section we discuss the processes affecting the formation 



and destruction of water molecules in disks. The discus- 
sion is structured as follow: first we explain how water 
molecules form in the disk atmosphere (Sec. I5.1|) : then 
we address the water (and OH) destruction processes 
and outline the main differences between T Tauri and 
Herbig AeBe stars (Sec. 15.21) . Finally we examine how 
non-stationary processes might affect the water content 
in different regions of the disk (Sec. 15.31 and 15 .4[) with a 
final discussion on the origin of the OH disk geometry 
correlation fSec. 15. 50 . 

5.1. Formation of water in disks 

The formation of water molecules in the warm disk 
atmosphere is a three steps process 

H + H + dust ->• H 2 + dust (6) 
+ H 2 -> OH + H (7) 
OH + H 2 ^H 2 + H (8) 

The formation of H 2 on warm dust grains in the 
disk atmosphere (Eq. El is justified by the findings of 
iCazaux fc Tielensl (|2002j) who measure a moderate H 2 
formation rate on warm (up to 900 K) dust grains. Fur- 
ther H 2 molecules might form through gas phase reac- 
tions fe.g. IGlasseold et al.ll2009h 

H- + H -> H 2 + e (9) 

The neutral-neutral reactions (Eqs. [7] & |5]) require high 
temperature to occur (T > 300 K). If H 2 is absent and 
hydrogen is mainly atomic, water may form through ra- 
diative association reactions (e.g. Kamp et al. submit- 
ted) such as 

H + 0->OH + /if (10) 
OH + H^H 2 + /ii/ (11) 

In both cases (neutral-neutral or radiative association 
reactions) the gas phase formation of water in the disk 
atmosphere strongly depends on the gas density. 

Recently a number of stationary disk chemical models 
have shown that water c an be efficiently f ormed in the 
warm disk surface layer. [Glassgold et al.l (|2009f ) (here- 
after G09) investigated in situ water formation with a 
model based on X-ray (only) heating and ionization of 
the disk atm osphere which is mo st appropriate for T 
Tauri stars. IWoitke et all (|2009fl (hereafter ProDiMo 
model) compute the thermo-chemical structure of Her- 
big AeBe disks wh ich generally have Lx/Lpuv << 1 
(jKamp et all [2008) using only UV/optical irradiation. 
According to their model there exists a hot water layer 
at distances between 1 - 30 AU and relative height z/r 
< 0.1 — 0.3 where water molecules are thermally de- 
coupled from the dust (T(H 2 0) > T(dust)). In this 
model OH (and CO) are abundant within 30 AU from 
the star and above an Ay of a few where water is not 
efficiently formed. Water is only found in deeper layers 
where it is shielded from photodissociation and densi- 
ties/temperatures are high enough to form it through 
neutral-neutral gas-phase chemistry. 

5.2. Destruction of water in disks 



Water depletion in Herbig AeBe disks 



9 



Pre-main-sequence stars emit strong ultraviolet radia- 
tion. For classical T Tauri stars UV photons are pro- 
duced by three main components: the stellar photo- 
sphere, the enhanced chromospheric activity, and the 
magnetospheric accretion. In an X-ray irradiated T 
Tauri disk (as the one investigated by G09) water 
molecules are dissociated in the disk atmosphere through 
charge transfer reaction with H + . In the case of early 
spectral type stars (B and A) as the ones studied here, 
the stellar photosphere is the major source of FUV radi- 
ation and overwhelms any contribution from mass accre- 
tion or stellar activity. The UV radiation field imping- 
ing onto the surface of the circumstellar disk affects the 
physical properties (such as the gas and dust tempera- 
ture) of the disk as well as its chemistry. In the case of 
the OH/H2O chemistry, the strength of the UV field reg- 
ulates the formation/destruction rate. Due to the strong 
soft UV radiation wate r can be easily photodissociated 
(e.g. IWoitke et al1l20"09l) 

H 2 + /ii/->OH + H (12) 
-^o + 2H (13) 
^0 + H 2 (14) 

Water can also be destroyed by charge transfer reac- 
tions as in the case of T Tauri stars (e.g. G09). However, 
given the intense ultraviolet radiation, UV photodissoci- 
ation is likely to be the major dissociation mechanism in 
the atmosphere of an Herbig Ae/Be star as we will show 
next. 

The minimum energy required to disso ciate water 
molecules is 5.1 eV (e.g. lHarich et al.l l2000h which cor- 
responds to a radiation of A = 2430 A. The photodisso- 
ciation cross section (k p d) of water varies with the color 
impinging radiation field. As an example, the cross sec- 
tion for an incoming radiation of a 10,000 K blackbody 
is nearly 4 times larger that for a 4,000 K blackbody 
(jvan Dishoeck et aLl 120081 ). This is due to the stronger 
soft UV radiation (relevant for the photodissociation of 
water molecules) of Herbig AeBe stars compared to T 
Tauri stars. Column 4 of Table [5] lists the UV luminos- 
ity of the program stars integrated between 1100 A and 
2430 A. The integrated luminosity is multiplied by a fac- 
tor of 10 - 4 '- 417 ™ to correct for interstellar extinction. 
The quantity An^o is the extinction at A = 1770 A and 
it is measu red from Ay u sing t he R = 3.1 extinction 
relation of ICardelli et al.l ([1989). The UV luminosity 
of the program stars ranges between ~ 1 — 10 3 L Q . In 
the case of a classical T Tauri star, the UV luminosity 
ranges between ~ 0. 1 — I Lq. In particular, for the two 
stars in ISalvk et al.l (pOOl (AS 205A, DR Tau) L uv is 
~ 0.1 — 1L . Thus, the UV radiation emitted by the 
Herbig AeBe stars in this sample is up to 4 orders of 
ma gnitude larg e r than that in the T Tauri stars studied 
by ISalvk et "aLl (|2008f) . Photodissociation is a plausible 
mechanim to explain the lack of hot water vapor lines in 
Herbig AeBe disks. 

A way to further test this scenario would be to de- 
tect OH ro-vibrational lines in the mid-infrared at 
high J rotational levels. In fact OH molecules formed 
by the photodissociation of water (eq. [TSl are yibra- 
tionally and rotationally excited (e.g. iDutuit et al.lll985l ; 



Ivan Harrevelt fc van Hemerti 120001: lHarich et al.l 120001) . 

As a consequence, the high rotational lev e ls of OH are 
easily populated (e.g. ICarringtonl Il964t lHarich et al.l 
12000 IBonev fc Mummal 120061) and mid-infrared spectra 
could be able to detect the high J- value transitions such 
as found in t he outflow of HH 21 1 (|Tappe et al.ll2009D and 
in TW Hya (jNaiita et al.ll2010D . A few high J OH lines 
have been detected in the mid-infrared towards Herbig 
AeBe stars with disks (Sturm, private communication). 
We note however that, depending on the density of the 
environment, the OH molecules might be thermalized 
very fast. This might prevent us from detecting any trace 
of water photodissociation. In this regard, the rotational 
diagram of L-band O H lines in the two stars studied by 
iMandell et al.l (|2008l) is characterized by a single rota- 
tional temperature. 

In addition to photodissociation there arc effects that 
could prevent us from detecting water vapor emission 
from the atmosphere of Herbig AeBe disks. The first ef- 
fect to consider is the temperature di fference between OH 
and H 2 lines detected here and in Sal vk et aLl ([2008). 
The lower energy state of the OH P4.5 doublet is ~ 
500 K while the energy state of H2O transitions covered 
by our observations are > 1000 K. Thus the OH emission 
potentially traces colder gas than the H2O emission. One 
might speculate whether the H2O non-detections are due 
to the different temperature layers probed in the disk. 
We find this possib ility unlikely given th at the OH tran- 
sitions detected by Mandcll et al.l (|2008D in AB Aur and 
MWC 758 have energies up to ~ 2360 K (OH P9.5), sim- 
ilar to the H 2 detected bv lSalvk et al.l (|2008l) . Thus the 
high OH/H2O line flux ratio in Herbig AeBe stars is a sig- 
nature of water depletion in the disk atmosphere rather 
than of transitions tracing different temperatures. 
Another argument that is often used to explain the 
paucity of molecular emission lines in Herbig AeBe disks 
is that the strong infrared excess might be able to veil 
the faint emission of mol ecular lines in the infrared (e.g. 
IPontoppidan et aLl I2010D . In this case the molecular 
emission lines are veiled by the dust continuum. In this 
regard we note that the we detect OH emission at sim- 
ilar wavelengths and similar sensitivity of the H2O ro- 
vibrational transitions detected in T Tauri disks. Thus, 
if water vapor is present, it must be located deeper in 
the disk (at higher Ay) where is colder than the OH de- 
tected here and/or thermally coupled with the dust (i.e. 
T(H 2 0) = T(dust)). 

5.2.1. OH photodissociation 

The most probable photodissociation channel of wa- 
ter is the one that produces an OH molecule (~ 80% - 
90% e.g. ICrovisieitfl989l: ICombi et al.ll200l lHarich et al.l 
12000) that is followed by 

OH + lw^O + H (15) 

The photodissociation cross section of H2O is only twice 
as large as that of the OH for an incoming black- 
body radiation of 1 0,000 K which is represen tative for 
an Herbig Ae star (|van Dishoeck et al.l 120081 ) . In ad- 
dition, the minimum energy needed to break the OH 
bound is 4.47 eV which corresponds to a radiation of 
2616A ( van Dishoeck fc Dalgarnol fl98l . very similar 
to that of water. This suggests that OH molecules 
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should also photo dissociate as is seen in comets (e.g. 
iCombi et a l. 2004). The dissociation of OH by UV 
photons fEq ll5|) is thought to be the main reservoir 
of the excited Oxygen ( 1 D) atoms which produce the 
6300 A transition often in solar system comets and in 



Herbig AeBe stars (e.g. [van Dishoeck fc Dalearnol 



Storzcr & Hollcnbach 1998; Morgcnth aler et al.1 



1984; 



2001 



Acke et al . 2005^3- The similar radial d istributions of 



OH and [O I] in HD 259431 (Sec. U and iBagnoli et all 
2010) might be a direct evidence of the formation of 
0( 1 D) atoms via OH photodissociation. 

The previous discussion suggests that the atmosphere 
of Herbig AeBe disks is depleted in water vapor. The 
detection of OH emission and the column density ratio 
OH/H2O > 1 suggests that OH molecules are produced 
again in situ (Eq. [7]). 

5.3. Transport of water in disks 

Non-stationary processes might affect the content of 
water in disks. Within a protoplanetary disk gas and 
solid particles c an migrate either inw ards or outwards 
(for a review see iCiesla fc C uzzi 2006j) or mixed by e.g. 
turbule nt motions. Outward s migr ation was suggested, 
e.g., by iStevenson fe Luninel (|1988|) to explain the for- 
mation of Jup it er. I nward migration was studied by 
ICiesla fe Cuzzil ([20061 ) with the aim at addressing the 
abundance of hot water vapor found in some T Tauri 
stars. Water can efficiently form in the disk interior on 
the surface of dust grains. If water-rich bodies such as 
icy planetesimals migrate inward and pass the snow line, 
the water ice evaporates from their surface and enrich the 
inner disk of water vapor. For such a process to happen 
however, the migration of planetesimals must be very fast 
and occur on a timescale shorter than the disk lifetime. 
Testing this hypothesis is important to probe the initial 
conditions of planetary systems. In this regards, two pa- 
rameters are important: 1) the abundance of d euterated 
water and 2) the ortho-para ratio of water (e.g. lEncrenad 
2008). Both these parameters are sensitive to the tem- 
perature formation and hence to the region where water 
forms in the disk (atmosphere or disk interior). We in- 
vite the intereste d reader to look into the recent work of 
iThi et al~l (|2010| ) for a detailed theoretical treatment of 
the formation of de uterated water. 

The case studied bv lCiesla fe Cu zzi (2006) is appropriate 
for T Tauri stars and does not include photodissociation. 
In the case of Herbig AeBe stars however, the transport 
and/or mixing of volatiles might induce the depletion 
of water in the disk interior by dredging up molecules 
to the surface layers where molecules are photodissoci- 
ated. Two timescales are important in this regard: 1) 
the chemical relaxation timescale {r c hem) which regulates 
how fast chemical reactions occur, hence the formation 
of water and 2) the mixing timescale. If the vertical mix- 
ing is faster than the chemical relaxation, the disk inte- 
rior can be depleted of water. Estimates of T c hem ior an 
Herbig Ae disk computed with ProDiMo are < 10 2 yr. 
The simulat ions of the transport processes in disks by 
lllgner et al.l (e.g. I2004D suggest that vertical mixing is 
slower and does not affect the global chemical evolution 

15 We note that another so urce of O^D) atoms is the direct 
photodissociation of water (Eq. I14j l 



of the disk. 

In this regard, ProDiMo predicts very different chemical 
structures between T Tauri and Herbig AeBe disks: T 
Tauri disks are colder and r c ^ em may be as long as 10 8 yr 
in the dense mid-plane where water is condensed in ice. 
Disks around Herbig AeBe stars instead have hotter and 
more (chemically) active environment and there is only 
very little amount of icy water in the mid-plane. 

5.4. Gas-dust decoupling and dust settling 

An important parameter for the photochemistry of 
disks is the physi cal decoupling of gas an d dust in the disk 
atmosphere (e.g. IMeijerink et al.ll2009| ). Direct observa- 
tional evidence of gas and dust physical de coupling was 
found for the Herbig Ae stars HD 101412 ([Fedele et al.1 
120011 and HD 95881 (IVerhoeff et al.ll2010h . In the case 
of physical decoupling, the disk atmosphere is depleted 
in large dust grains (which settle into the midplane) and 
water molecules are unshielded from the UV radiation. 
In such a case the molecular content in the disk atmo- 
sphere might be easily reduced by UV photodissociation. 
The situation changes dramatically if water and OH self 
shield ing is taken into account as in iBethell fe BergirJ 
(2009). According to their model, in a dust-dcpletcd 
disk atmosphere water and OH molecules become the 
major source of UV opacity and are able to block the 
photodissociative radiation from penetrating further into 
the disk. As a consequence, the abundance of water is 
enhanced even in the presence of a moderate FUV ra- 
diation (Lpuv = 10L (T)). One problem of th e water 
self-shielding model of IBethell fe Berginl (|2009| ) is that 
it assumes the presence of H2 in the dust-depleted disk 
atmosphere. Without dust however, H2 form via the 
H route which produces an overall lower abundance of 
molecular hydrogen. On the other hands, the model of 
IMeijerink et al.l (|2009fl as well as G09 and ProDiMo do 
not include water self-shielding. This makes difficult to 
compare the observations with predictions from different 
models at this stage. 

We do not have direct evidence of physical decoupling 
in the systems studied here but we do have evidence for 
thermal decoupling. If we assume a gas tem perature of 
T - 700 K as found bv lMandell et al.1 ([20081 ) for AB Aur 
and MWC 7 58, the gas is hotter than t he dust (T(dust) 
< 300 K, e.g. lKamp fc Dullemondll2004T ) at a spatial scale 
of ~1 - 20 AU, the size of the OH emitting region (Sec. 
[4]). Thus, the OH vapor detected here is likely thermally 
decoupled from the dust. 

5.5. OH emission and disk geometry 



In Sec. 14.11 we found that the presence of OH ro- 
vibrational lines correlates with the geometry of the disk: 
OH emission is mainly detected towards flaring disks 
(group I). As we saw in Sec. 13.1.11 the excitation of the 
OH ro-vibrational transitions depends on the tempera- 
ture and density of the gas in the case of thermal emis- 
sion and on the intensity of the UV and near-infrared 
radiation field in the case of fluorescence. In the case 
of Herbig AeBe stars the disk temperature is regulated 
by the UV output of the star as well as the near-infrared 
continuum emitted by the disk. Thus the intensity of OH 
ro-vibrational lines depends on the UV radiation field im- 
pinging into the disk surface regardless of the excitation 



Water depletion in Herbig AeBe disks 



11 



mechanism. In a flaring disk geometry, the disk surface 
area that is illuminated by the central star is much larger 
than in a self-shadowed geometry. In the latter case OH 
might still be excited but the emitting area is small (and 
the emission weak) and confined to the inner part of the 
disk. 

A trend is also found between CO fundemantal ro- 
vibrational emission at 4.7/im and disk flaring (van der 
Plas et al. in preparation). CO emission originates from 
a larger disk area in flaring disks compared to flat disks. 
A larger sample of Herbig AeBe disks and more OH tran- 
sitions are necessary to pin down the OH line excitation 
mechanism. 

6. CONCLUSIONS 

This paper presents high resolution L-band spectro- 
scopic observations of Herbig AeBe stars with disks. 
Among the 11 stars analyzed here, four show hot OH 
emission which appears to correlate with the disk geom- 
etry and, in turn, with the UV irradiation impinging on 
the disk. The detection rate of OH emission is 70% for 
disks with a flared geometry (in cluding the two systems 
analyzed bv lMandell et alJ l2008) and - 13% (1 out of 6) 
for self-shadowed disks. The OH ro-vibrational lines are 
spectrally resolved allowing us to measure an extension 
of ~ 10-30 AU for the OH emitting region. In contrast 
to T Tauri, Herbig AeBe stars do not show evidence of 
hot water vapor in their spectra. The detections of OH 
emission lines and the observed OH/H 2 column density 
ratio (> 1) indicate that the atmosphere of disks around 
Herbig AeBe stars is depleted in water molecules. Given 
the stronger UV radiation field of Herbig AeBe stars com- 
pared to T Tauri stars, a plausible explanation for the 
non-detection of water lines is that water in the disk at- 
mosphere is dissociated by UV photons. The absence of 
hot water vapor in the disk atmosphere does not preclude 
the existence of colder water deeper in the disk (at lower 
temperature and higher Ay). Far-infrared water lines 
have been detected with the Herschel Sp ace Observa- 
tory t owards the Herbig Ae star HD 100546 (jSturm et al.1 
2010) (although further confirmation is needed). This in- 
dicates the presence of warm water vapor in the disk. 
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APPENDIX 

ON THE ASYMMETRIC OH LINE PROFILE IN V380 ORI 

As we saw in Sec. I4.2[ the OH line profile in V380 Ori is asymmetric, with the blue-shifted component brighter 
than the red-shifted one. A similar asymmetry was found towards other protoplanetary disks in the [OI] 6300 A line 
for the H erbig Ae star HP 100546 (|Acke fc van den~A nckcr 20061) and in the CO fundamental rovibrational lines in 
EX Lupi ()Goto et al.ll201ll) . Such an asymmetry is either due to a deviation from the Keplerian motion of the gas 
or to a non homogeneous distribution of the emitting gas. The inner radius of t he OH emitting region (~ 2AU) is 
too large to be coincident with the dust sublimation radius. Alecian et al. (2009) found evidence of a close low-mass 
star companion to V380 Ori. They determine a projected angular separation asin(i) < 0.33 AU. Regardless of the 
inclination, the orbit of the companion is certainly within the inner rim of the OH emitting region. Circumbinary 
disks have inner gap of the order of ~ 2 - 3 times the binary separation (|Artvmowicz fc Lubowlll994D which brings 
the truncation radius to ~ 1 - 2 AU for V380 Ori (depending on the eccentricity, inclination and mass ratio of the 
binary). This is remarkably close to the OH inner radius estimated here. The companio n might be also the source 
of the perturbation of the gas dynamics seen in the OH line profile. Recent simulation by I Re gal v et al.l (|2011l ) show 
that in the case of a circumbinary disk the velocity distribution of the gas differs from the circular Keplerian case. 
The disk may become eccentric (due to tidal interaction with the lower mass companion) and the velocity profile of 
the emerging lines is asymmetric in a fashion similar to Fig. [2j 
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Table 4 

OH and H2O lines parameters 



Star OH FWHM OH |EW| L OH Lh,0 @ 2.931 /jm log(N(OH)) log(N col (OH)) log(N(H 2 0)) log(N co i(H 2 0)) 

[km s _1 ] [10 _6 /im] [10 -15 erg cm -2 s -1 ] [10 5 erg cm -2 s -1 ] [molecules] [mol./cm -2 ] [molecules] [mol./cm — 2 ] 



UX Ori 








< 3.0 




< 44.3 


< 15.4 


< 44.3 


< 15.5 


HD 34282 








< 2.0 




< 43.2 


< 14.3 


< 43.2 


< 14.4 


CO Ori 








< 1.0 




< 43.8 


< 15.0 


< 43.9 


< 15.0 


V380 Ori 


27 


8 


8.4 (± 5.0) 




< 4.5 


45.0 


16.1 


< 44.5 


< 15.6 


BF Ori 






< 1.0 




< 43.8 


< 15.0 


< 43.9 


< 15.0 


HD 250550 


18 


2.8 


1.4 (± 0.1) 




< 2.0 


44.6 


15.9 


< 44.5 


< 15.8 


HD 45677 






< 10 




< 43.9 


< 15.0 


< 43.9 


< 15.1 


HD 259431 


45** 


12 


14.5 (± 2) 




< 1.0 


45.7 


15.9 


< 44.3 


< 14.5 


HD 76534 






< 1.5 




< 44.5 


< 15.7 


< 44.6 


< 15.7 


HD 85567 


30** 


8 


10.0 (± 30) 




< 3.0 


46.1 


17.3 


< 45.4 


< 16.5 


HD 98922 






< 1.2 




< 44.6 


< 15.7 


< 44.6 


< 15.8 



Upper limits arc computed as the product 3 X tr x AA, with a standard deviation between 2.930 fj,m and 2.936 /im, and AA = 30 km s -1 . The 
OH line flux refers to the average flux of the two transitions of the 2 n 3/2 P4.5 doublet. The FWHM and EW are the moan of the FWHM 

and EW of the two transitions. Total and column densities arc computed assuming optically thin emission, LTE and T(OH) — T(H20) — 700 K 
and T(OH) = T(H 2 0) = 1000 K respectively. 

In this case the OH doublet is blended (Fig. \5\ & 0. The FWHM and EW arc computed as the half of the total doublet FWHM and EW. 



